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Abstract 

Groups of animals often tend to arrange themselves in flocks that have char- 
acteristic spatial attributes and temporal dynamics. Using a dynamic contin- 
uum model for a flock of individuals, we find equilibria of finite spatial extent 
where the density goes continuously to zero at a well-defined flock edge, and 
we discuss conditions on the model that allow for such solutions. We also 
demonstrate conditions under which, as the flock size increases, the interior 
density in our equilibria tends to an approximately uniform value. Moti- 
vated by observations of starling flocks that are relatively thin in a direction 
transverse to the direction of flight, we investigate the stability of infinite, 
planar-sheet flock equilibria. We find that long-wavelength perturbations 
along the sheet are unstable for the class of models that we investigate. 
This has the conjectured consequence that sheet-like flocks of arbitrarily 
large transverse extent relative to their thickness do not occur. However, 
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we also show that our model admits approximately sheet-like, 'pancake- 
shaped', three-dimensional ellipsoidal equilibria with definite aspect ratios 
(transverse length-scale to flock thickness) determined by anisotropic percep- 
tual/response characteristics of the flocking individuals, and we argue that 
these pancake-like equilibria are stable to the previously mentioned sheet 
instability. 

Keywords: Flocking, Swarming, Stability Analysis, Biological Aggregation, 
Density Profiles, Continuum Flock 
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1. Introduction 

The formation and movement of groups of animals is a collective phe- 
nomenon emerging from behaviors of the individual group members. One 
particularly striking example is the coordinated motion of large flocks of 
birds (such as starlings The goal of this paper is to utilize a simple 

model to study equilibria and stability of flocks. In particular, we note the 



following quote from Ref. 
ling flocks: 



1] which reports an empirical investigation of star- 



"Perhaps the most interesting morphological result is that flocks 
seemed to have a characteristic shape, being thin in the direction 
of gravity and more extended perpendicular to it. The propor- 
tions of the flock were well defined, with only weakly fluctuating 



2 



aspect ratios, despite showing a wide range of sizes. Our abil- 
ity to conclude this stems entirely from the fact that we were 
able to analyze several groups with very different sizes (dimen- 
sions and number of birds). Nonspherical shapes have also been 
observed in fish schools, the average proportions Ii.I^-Iz ranging 
from 1.0:1.7:2.1 in pilchards, ...to 1.0:3.0:6.0 in saithes, ...and 
1:3:4 in herrings, .... These values are comparable to those we 
found for starlings: 1.0:2.8:5.6. 

Analysis, to be presented here, shows that our simple model reproduce these 
features. In particular, a flock tends to have the same aspect ratios indepen- 
dent of flock size. 

There are two main modeling paradigms that have been used to inves- 
tigate collective behavior of flocks. In one, the position of each flock mem- 
ber and its interactions with other flock members and the environment are 
evolved 2j-|9|, while the other is based on a continuum model evolving the 
density of individuals (g[ [j-ol— Tl7 1 . In addition, within the class of continuum 
models, two subclasses may be distinguished which we call kinematic and 
dynamic (after Ref. jlO]). In kinematic models the macroscopic flock veloc- 
ity v(x, t) is instantaneously determined by the flock density, while dynamic 
models employ a velocity evolution equation specifying dv/dt. In what fol- 
lows, we will use a dynamic continuum model. A main motivation for our 
use of a continuum description is that continuum models are better suited 
to analyze stability and wave propagation. For a review of this and other 
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collective behavior phenomena, see Refs. 18M23 1 . 

This paper is organized as follows. In Sec. |2] we present our continuum 
flock model. In Sec. [3] we apply our continuum model to obtain planar, one- 
dimensional equilibrium solutions, as well as ellipsoidal three-dimensional 
equilibrium solutions. We show that certain conditions on the model equa- 
tions lead to finite flocks (i.e., zero density of flock members outside a finite 
spatial region) and others do not. Similarly, we find that, if the model is fur- 
ther constrained, the flock's interior density does not blow up with increasing 
flock size (an unrealistic situation), but rather tends towards a well-defined, 
approximately uniform interior density. In addition, we discuss how our re- 







24 


i 



the 'topological distance' as advocated in Ref. |24j) 1 l. In Sec. H]we investigate 
an instability of a one-dimensional ('sheet') equilibrium to long-wavelength 
perturbations along the sheet. Conclusions and discussion are given in Sec. 
5. 



1 By topological distance, Ref. 



24j means that a flock member interacts with a fixed 



number of nearest flock members, independent of their geometrical distance. 
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2. Model Formulation 



2.1. Continuum Model 



The equations for our model are similar to those of Ref. 25], 



dp 
dt 



+ V ■ (pv) = 0, 



(1) 




VP(p) - W + - ( 1 -\ ) v 



(2) 



In Eq. (JTJ) , p(x, t) denotes the density of flock members at spatial point 
x = (x,y,z). Equation ([T]) represents the conservation of the number of 
flock members. The density p(x, t) may be thought of as the number of flock 
members in a small region of dimension / centered at the point x divided 
by the volume of the region. An inherent assumption is that the flock is 
large enough that regions of size I can contain many flock members yet still 
satisfy the condition that I is substantially smaller than the macroscopic 
characteristic spatial scales of the flock. The quantity v(x, t) denotes the 
macroscopic flock velocity at point x given by the average over a region 
of size / centered at x of the individual velocities of flock members. The 
quantity -P(p) will be referred to as the 'pressure' due to its analogy with the 
fluid dynamical pressure and the similarity of the form of Eq. (T5]) with the 
momentum equation of fluid dynamics. The term P(p) may be heuristically 
thought of as arising from dispersion, relative to v(x, t), of the velocities 
of individuals within the regions of size I, referred to above, coupled with 
short-range repulsive behavior of individuals, where this short-range behavior 
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is motivated by collision avoidance! In what follows, we shall investigate 
the consequences of assuming different forms for the function, P(p). The 
quantity C/(x, t) is a potential modeling the longer-range attractive behavior 
of individuals necessary to form a flock. The last term in Eq. represents 
the tendency of flock members to have a preferred speed vq with respect 
to the medium through which the flock moves (e.g., air and water in the 
cases of flying birds and swimming fish), and r is the rate of relaxation of 

nnn 

v(x,t) = |v(x, t)\ to this preferred speed (see, for instance, Refs. [13l. |26|-|28|). 
Here we choose U to be 

U = J u(x - x')p(x') d D x', (3) 

where d D x' represents a .D-dimensional differential volume element (D = 1 
or 3 in this paper). For convenience, the function u(x) will be taken to have 
the particular form given by the solution of the equation, 

V^m(x) - k 2 m(x) = u 5(x), (4) 

where, 

= V ■ K ■ V (5) 



^.2 ^.2 

IK = SiSi— + s 2 s 2 — + s 3 s 3 — , (6) 

Hil Hie) K>n 



2 As in Ref. [25| an additional term, representing the tendency of flocking individuals 
to align their directions of motion, v/|v|, with those of their neighbors, can be included 
on the right-hand size of Eq. ((2J. However, this term turns out to be identically zero for 
the equilibria we investigate, and, for simplicity, we have thus omitted it. 
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and k, ^1,2,3, and Uq are constants, and {si, s 2 , s 3 } are a mutually orthogonal 
set of unit vectors. Note that V 2 is an anisotropic form of the Laplacian 
operator; V 2 = V 2 for k 2 = n\ = k 2 = ft 2 . This form for V 2 is introduced to 
allow for the natural anisotropy of behavior and visual perception of flocking 
individuals. For example, in flying birds, Si(x, t) = v/v, s 2 is the average 
direction of orientation from wingtip to wingtip within regions of size I, and 
s 3 (x, t) is the average orientation of the direction perpendicular to the plane 
formed by the head and outstretched wings (i.e., s 3 is typically in the di- 
rection of gravity). For example, the solution of Eq. (j3J) for u(x) in three 
dimensions (D = 3) is 

1 /2 

u(x) = [u /(47rr)] exp(— nr), r = jx ■ KT 1 • x} , (7) 

where we have rejected the solution to Eq. that blows up at r — > oo. 
Here, k~ 1 represents a chosen reference length scale for interactions between 
flock members. 

We note that, while we have accounted for anisotropy in flock mem- 
ber behavior and perception, we have not attempted to account for for- 
ward/backward asymmetry. In particular, one expects an individual to be 
more sensitive to conditions in front of it than to conditions at an equal 
distance behind it. This type of forward/backward asymmetry is often ac- 
counted for in models that are based on equations of motion for each of many 
simulated individual flock members (for instance see Refs. jj], Q]); but, so 
far, has not been incorporated in a continuum model. 
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2.2. Other Continuum Models 



Continuum models have 



)een used for many 



behavior of animal groups [8 



' [io-12, Q. 



15 



17 



years to mode 



25 



27 



28 



30 



collective 



435|] . Most 



such models use continuity of agents to constrain the motion of individuals. 

HQ, 



In addition to this, some researchers 



17 



34| set the velocity to be a 



particular function of the density and velocity and their spatia 



These models, which could be called kinematic (after Ref. 



12 



25 



28 



30 



31 



derivatives. 



10(| ) , a re to be 



35| such as 



distinguished from dynamic models 

an 

the model just presented. In dynamic models the macroscopic flock velocity 
evolves in time through an equation for dv/dt (as in our momentum-like 
equation given in Eq. <^). 

Similar to discrete models, continuum models employ a subset of terms 
that represent the attraction, repulsion, orientation, self-propulsion, and 



ed the anisotropic behavior 

2m. 



noise. Only a few continuum models have mode 
of sensing and response, most notably Refs. [12I . 

Repulsion of nearby individuals is sometimes considered a nonlocal term 
[ill . Q, analogous to the attractive potential U given in Eq. ([31), but with 
the opposite sign and a shorter range. Since the biologically relevant limit 
is for repulsion to have a shorter length-scale, we have directly incorporated 
repulsion through the local gradient of a 'pressure' in Eq. ([2]); this is similar 



to Ref. 23|. Attraction is longer range and is included in our model in a 



nonlocal manner (Eqs. ([3]) and (jl]) for U and it), and nonlocal modeling of 
attraction is also done in most other continuum and discrete flocking mod- 



els. Some authors have made fully local models (like Refs. 



31 



32|), 



but these models cannot accurately simulate long-range attractive behavior. 
Some researchers Q, [ill, 12, 25, 28, 34] use a nonlocal convolution of an ex- 



ponential kernel and the density, similar to our model. Others use attractive 
terms with kernels including piecewise continuous and power-law dependen- 



cies (see, for instance, Refs. 



13, 



34 



36|). To what extent the form of these 



attraction terms affect the model and conclusions is not clear. 



3. Equilibria 

We obtain solutions to our model equations that represent a uniformly 
translating flock. To this end we set 

v = v s, (8) 

where vq and the unit vector s are constant everywhere within the flock. 
After translating to a frame moving with the constant speed v, Eqs. ([2]), ©, 
and (jlj) yield the equilibrium equations, 

--VP(p)-W = 0, (9) 
V 2 p U - k 2 U = u oP , (10) 

where Eq. ( fit)]) results from Eq. (J3]) by replacing x by x — x', multiplying by 
p(x') and then integrating the result over x'. Note that Eq. ([9]) applies only 
inside the flock (p(x) > 0), since Eq. fl2]) is irrelevant for p = (outside a 
flock). Equation (Q can be rewritten 

V [f(p) + U] = 0, (11) 
9 



where 

f( P )= r^dP(p'). (12) 

Jo P 

Thus within the flock 

C/(x) + f(p) = U B , (13) 

where U-q is a constant. Since we seek equilibria with p = on the flock 
boundary, we impose / = and U = Ub on the boundary. In principle we 
can solve Eq. (TT3~]) for p, yielding a result of the form 

p = F{U B -U), (14) 

where 

= x. (15) 

Combining Eqs. (11 4p and (TTOT) yields a nonlinear modified Helmholtz equation 
for £7, 

f F(Ub - U) inside the flock, 
Vjf/-K 2 t/ = Mo < (16) 
I outside the flock. 

Integration of Eq. f[T5"|) across the flock boundary yields the condition that U 
and VU are continuous across the flock boundary. Outside the flock (p = 0) 
the magnitude of the relevant solution of Eq. (fl6l) decays monotonically to 
zero at infinity in accordance with Eqs. (JTj) and ([3]). Also, the potential must 
be confining within the flock to balance the repulsive action of the pressure. 
These facts imply that U is negative everywhere. 

We will consider a planar ('sheet') and three-dimensional ellipsoidal ('pan- 
cake') equilibria. In the case of 'sheet' equilibria, we assume a one- dimensional 
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density profile, and thus we refer to this as a planar one-dimensional solution. 
In this case U depends only on x and Eq. (1TB1) reduces to 



,2#U 2tt F{U b -U) inside the flock, 

ax I outside the flock, 

where £ 2 = x-K-x = X)n=i( s n ' X ) 2 ( K / K n) 2 - in the case of three-dimensional 
ellipsoidal equilibria, without loss of generality, we can take Si = x, s 2 = y, 
and s 3 = z. We further assume that U is a function of the single variable 

r = {x • K" 1 • x} 1/2 = {(k iX ) 2 + (K 2 yf + (k 3 z) 2 } 1/2 /«. (18) 

In this case Eq. f lT6|) reduces to 

1 d ndUir) 2rr/ . f F(U B -U) inside the flock, 

— — r 2 — - K 2 f/ (r) = m t (19) 

r dr dr I outside the flock. 

3.1. Finite Flock Condition 

We now consider the effects of different choices for the 'equation of state', 
P = P(p). A main conclusion will be that if, as p — > 0, the pressure P(p) 
behaves like 

P(p) = cp\ (20) 
then finite equilibria are possible only if 

7 > I- (21) 

Here by 'finite equilibria' we mean equilibria in which p > is continuous in 
a region 71, going to p = on the region boundary d7Z, and p = outside 
1Z. Equation fl2Tj) follows as a consequence of the fact that Eq. f|T2|) implies 
f(p) ~ p 7_1 near p = 0, and thus / — > as p — > + only if Eq. (|2"T|) holds. 
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3.2. Behavior at the Flock Edge 

Let r]± denote a coordinate measuring the perpendicular distance to the 
flock edge, where r]± > corresponds to points inside the flock. Expanding 
Eq. OH]) about a point p on the flock boundary and assuming Eq. (I2(jp for 
small p, Eq. fl9]) yields 



1 dp^ 



pdrj A 



J, J 



1 dU 



c drjj 



and the dependence of p near the edge (r]± small) is 



i 

7-1 

P ~ vl 



(22) 



(23) 



as illustrated schematically in Fig. [TJ 






Figure 1: Schematic of the density dependence near the edge of the flock for different 7. 



3.3. The case, P(p) = cp 2 

In this case, Eq. f lT2|) gives f(p) = 2cp, which, when combined with Eqs. 
flU and dU, yields 

F(U B -U) = [U B - C/(x)] /2c. (24) 
12 



Thus, for this special pressure dependence, Eq. (TIB]) is linear in U(x), and 
this facilitates an analytical solution of the problem. In the one-dimensional 
planar case, the solution of Eq. (fTTj) within the flock (\x\ < xb) is 

p(x) cos(y/K — 1 k^x) — cos(y/K — 1 K^Xb) , < 

p(0y~ 1 -cos(v / F^I<" 1 xb) ' 

with = uo/(2ck 2 ). Outside the flock (|a;| > xb and p = 0), the solution for 

Eq. (JT7j) is a negative exponential, U(x) = Ub exp(— — Xb|). Use of 
the conditions at the boundary (that U is continuous and has a continuous 
derivative) determines the unknown xb in terms of the physical parameters 
uq and c, 



sin(VK - 1 KfVg) = - ViiT- 1 cos(VAT- 1< _1 ^b)- (26) 
Several solutions of this type are shown in Fig. 12(a). 



In the three-dimensional ellipsoidal case, the solution to Eq. fll9p is 

p(r) 



P(0) 



ktb sin (y/K - 1 rer) 
«r sin [Va — 1 ktb) 



1 



y/K — 1 /?r B 
sin (V-K" — 1 k? "b) 



- 1 



(27) 



Outside (r > tb), the density is zero and the boundary conditions give the 
boundary position as 

K 



sm(\/K — 1 m"b) = — V — 1 cos(v / K — 1 ktb) 



l-K- 



• (28) 



Solutions of this type are shown in Fig. [2](b) for various values of K. 

Note that, as a consequence of the linearity of Eqs. ifPTj) and (fl9i) . for 
7 = 2, the density profiles of p(x)/p(0) and p(r)/p(0) and the normalized 
boundary location, /t£ _1 a;B and ktb, are determined by the model parameter 
f^, and are independent of p(0) (e.g., this is not so for P(p) = cp 1 with 
7^2). 
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3.4- One- Dimensional Equilibria 

For a general pressure, P(p), the solution in one planar dimension can 
be formally obtained by quadrature. Multiplying Eq. ( ITTI) by dll/dx, Eq. 
( ITT)) can be written as an exact differential. Integrating this once yields 
a separable, first-order differential equation which, using Eq. (Jill), can be 




b) 



P/P(0) 




AT 



Figure 2: Density profiles for planar (a) and ellipsoidal (b) flock equilibria with P = cp 2 
and K — u /2ck 2 . 
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integrated to give 

x B -x = i / =, 29 

Ju(x) y/n 2 U 2 " 2u P [F(Ub - U)\ 

where 

x B = £ =, (30) 

A/ min v 7 ^ 2 - 2 «oP [^(% - U)\ 

and f/ min = U (0) is the root of 

K 2 Ul^ - 2u P [F(U B - U min )} = 0, (31) 

such that the integrand is positive between f/ min and Here is the 
value of the potential at the boundary of the flock and may be considered to 
be determined by specifying the number N of individuals per unit transverse 
area in the flock, 

/+oo 
p{x) dx. (32) 
-00 

Since the potential U is attractive, U and U-q are negative, and U-q is the 
maximum value of U. 

We now consider pressures of the form P(p) = cp 1 with 7 ^ 2. For 7 < 1, 
we can show that there are no solutions where p(x) > and p(x)dx < 00, 
and we therefore only consider 7 > 1. Using Eqs. (f!2l) and (flBl we obtain 

"7-I 



F(U B -U) = p 



cry 



(U B - U) 



7-1 

(33) 



Figure [3] shows density profiles for 7 = 1.5, 3, and 5.0 for P(p) = cp 1 . In 
Fig. El 

<*> 

15 




Figure 3: ID Equilibrium density profiles for different pressures. Density profile for a 
pressure of P = cp 1 for a) 7 = 1 .5, b) 7 = 3.0, and c) 7 = 5.0. Each profile has the same 
value for nN/{p£) w 54. 
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These three plots all have the same value of nN/(p£). Note that the results 
for 7 = 3 and 7 = 5 are qualitatively similar and show that p(x) rises 
rapidly as one moves into the flock from its boundary (x = xb), becoming 
approximately constant in the interior of the flock. In contrast, the result 
for 7 = 1.5 is quite different, being peaked at x — 0, with a very much larger 
value of p(0) and a very much smaller value of xb- 

*^Vmu./(P?) 

10 s 
10 4 
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Figure 4: This figure shows the minimum number of individuals that will support a finite 
flock for the case of a pressure P(p) — cp 1 for f < 7 < 2. 



When 1 < 7 < 2, there is a minimum number of individuals, parametrized 
by K p N/ (p£), that will support a finite flock. Fixing Ub sets the total number 
of individuals in the flock. A large \Ub\ corresponds to a large flock size. The 
smallest \Ub\ can be is 0, which corresponds to a minimum flock size. This 
minimum value of nN/{p^) is plotted as a function of 7 in Fig. HI We see 
from the figure that the minimum flock number approaches infinity as we 
approach 7 = 2 from the left. 
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In Fig. [5J we show density profiles for different values of nN/{p^) in the 
case 7 = 3. We see that, for small kN/ (p£), the density profile is not constant 
near the center of the flock (x — 0). However, as the flock number increases, 
we approach the flat internal density and sharp cutoff depicted in Fig. E^b). 

In Fig. [H] we plot the central flock density p(0)/p versus K,N/(p£) and the 
flock boundary k^^x-q versus nN/(p£) for 7 = 1.5 [Figs. E[a,c)] and 7 = 3.0 
[Figs. E](b,d)]. For 7 = 1.5, we see in Fig. Eta) that the density at x — 0, 
p(0), increases faster than linear as KiV/(p£) increases. In the case 7 = 2 
(described in Sec. 13.31) . we found that p(x)/p(0) and xb are independent of 
p(0), implying that p(0) is linearly proportional to N. The faster-than-linear 
increase of p(0) with nN/(p£) is consistent with Fig. [6](c) which shows that 
the position of the flock boundary, xb, decreases with KiV/(p£). For the 
case 7 = 3.0, we see in Fig. [6](b) that the density at the center of the flock 
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Figure 5: Density profiles for the pressure P(p) = cp 3 for various kN/(p£). The curves 
represent kN /(/?£) values of 4.8, 16.8, 36.7, and the last, 54, is the same curve as Fig.[3]Jb). 
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saturates at a constant value as nN/{p^) increases. Consistent with this, 
Fig. [6](d) shows the flock boundary position x-q increasing proportionally to 
kN/(p£) for large kN/(p£). Hence we conclude that for pressures of the 
form P(p) = cp 1 , with 7 = 3, large flocks of different sizes will have the 
same flat interior density. In fact we find this to be the case for any 7 > 2. 
Furthermore, for 7 > 2 we can obtain an explicit expression for the central 
density, as follows. From Eq. ( TT7l) . we see that, if the density is nearly 




kN/(p& 



Figure 6: Maximum density, p(0)/p, versus flock number, nN/(p£), for a) 7 = f.5 and b) 
7 = 3.0. Boundary position, k£,~ 1 xb, versus flock number, kN/(p~£), for c) 7 = f.5 and d) 
7 = 3.0. 
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constant in the interior of the flock, then d 2 

p( ) -—U n 



U/dx 2 = and 



(35) 



From Eq. (j5T|) 



« 2 ^. n _ 2mo c P (0)t = 0. 



(36) 



Using Eq. (|35|) in Eq. (I3TI) gives 



P(0) 



(37) 



Equation (l3Tj) is plotted as the solid curve in Fig. [7J the open circles in Fig. 
[7] are the central density for a large flock evaluated by numerical integration 
of Eq. 029]). 

As noted in the quote in Sec. (U it is common for animals of a given species 
(e.g., starlings) to form flocks with different numbers of flock members. As 
the number of flock members increases the spatial extent of the flock also 
increases. For a model with P(p) = cp' 7 , our results above show that this 
type of observed behavior is only consistent with 7 > 2. Furthermore, it 
is expected that the interior density of larger and larger flocks eventually 
saturates, since individuals would be expected to become 'uncomfortable' 
beyond a certain crowding density. This is also consistent with 7 > 2 since 
we obtain an interior density for large flocks, Eq. ( 1371) . that saturates with 
large kN/(p£). 

3.5. Three- Dimensional Ellipsoidal Equilibria 

Unlike our analysis of the planar one-dimensional case in Sec. I3.4[ the 
three-dimensional ellipsoidal problem, Eq. ( Fl9l) . cannot be solved by the 
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method of quadrature. While one could solve Eq. ( I19p numerically, we feel 
that this is not necessary, as insight from the planar case can be used to give 
a good understanding of the ellipsoidal case. In particular, we focus on the 
physically interesting case of 7 > 2 and large flocks. From our discussion in 
Sec. 13.41 we expect that the interior density will again be flat, and, neglecting 
r~ 2 d/dr(r 2 dll/dr) in the flat region, we again obtain Eq. (|37|) for the inte- 
rior density. Furthermore, consider the thickness Ar of the transition region, 
t-q — Ar < r < r^, where p(r) rises from p(r^) = (at the flock edge) to 
p(r) near p(0). For large flocks, we expect that Ar <C r-Q. Thus the effect 
of curvature of the flock boundary will have little effect on p(r) within the 
transition region, and the spatial dependence of p on r from the boundary 
will be nearly the same as in the one- dimensional planar case. Thus, for 

p(0)/p 
1.0 ■ 

0.8 ■ 

0.6 ■ 

0.4 ■ 

0.2 ■ 

123456789 10 

Figure 7: Plot of the density at x = versus 7 with a pressure of the form P(p) — cp 1 . 
The black curve is the theoretical expression Eq. (|37]). The open circles are the central 
densities of large flocks obtained by numerical integration of Eq. (|29|) . 
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large flocks p(r) will be the same as p(x) given Figs. |3](b and c) provided 
that we replace the horizontal axis variable n^~ l x in Figs. |3^b and c) by Kr, 
and provided that the planar and ellipsoidal flock sizes are the same (in the 
sense that x-q = re)- Finally, replacing r by Eq. ( ITS"]) we obtain the pancake- 
shaped ellipsoidal flock density as a function of (x,y,z). Referring to the 
quote in Sec. HJ the aspect ratios of the flock shape are l:(/t 2 /Ki):(/t3/Ki), 
and, in agreement with the observations, these aspect ratios for our model 
are independent of the flock size. 

3.6. 'Topological Distance' 
A recent paper 



24J | suggests that the interactions between flocking star- 
lings is through what they call the 'topological distance', rather than the 
geometric distance. Thus, for example, in the case of isotropic interactions, 
each bird may be thought of as interacting with a fixed number of other 
birds that are closest to it independent of the geometrical distance between 
them. In contrast, our analysis above has assumed that interaction strength 
between individuals falls off with increasing geometrical distance. In this 
subsection we show how our analysis above may be simply adapted to the 
case where flocking individuals interact through topological distance. We 
accomplish this by replacing the operator V^, defined by Eq. fl5]), by the 
modified form 

= ■ (k^\ ■ V, (38) 

where p is the flock density, and is a somewhat arbitrary reference density 
introduced to leave the units of as the inverse of length squared. A 
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convenient choice for might be the density in the center of the flock. We 
now consider the effect of this modification in the case of a one-dimensional 
planar flock for which our formulation using geometrical distance yields Eq. 
(IT7|) . With our modification, Eq. (|T7|) still applies, but with x replaced by 

x(x) = P ^ dx', (39) 
Jo P* 



which is essentially the topological distance, as described in Ref. [24j. Consid- 
ering a large flock, the central density p(0) is thus still given by our previous 
analysis, Eq. ( 137|) . and the density for large flocks is again approximately 
constant in the flock interior. To find the dependence of p on distance from 
the flock edge we replace r)± in Eq. fl23|) by the corresponding topological 
quantity, fj±. This yields 

P-V 1 /^, (40) 
which is to be contrasted with the dependence p ~ r} 1 /^ 1} in Eq. ([23]). 



Note that in Sec. 13.41 we have argued that 7 > 2 for physical flocks with an 
assumed pressure of the form P ~ p 1 . Thus, in this case, Eq. f T4"0"j) continues 
to predict that p goes continuously to zero at the flock edge. Furthermore, 



by the discussion in Sec. 13. 5[ our consideration above for the one- dimensional 
planar case, directly translates to the case of an ellipsoidal flock. 

3.7. Other Continuous Equilibria 

In our approach we regard density profiles that go continuously to zero at 
the flock boundary and approach a constant uniform interior density at large 
flock size as 'reasonable', and we have delineated conditions in our model 
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class where this holds. Furthermore, our work shows that our reasonable 
equilibria are robust in the sense that they exist in regions of parameter 
space (rather than just at a critical parameter value). 

In contrast, Ref. presents equilibria where the density profile goes 
discontinuously to zero at the flock edge, while another [8i] finds equilibria 
that give continuous density distributions that asymptote to a small constant 
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34j obtain profiles of 



density far from a localized flock. References [10, 
the type that we regard as reasonable; they, however, do so using a kinematic 
model rather than a dynamic model as used here (see Sec. 12.21 for our defini- 
tion of dynamic vs. kinematic). Our flock solutions, in which the density is 
strictly zero outside the flock, is similar to 'compacton' solutions of certain 
nonlinear PDEs [37J. 

We note that we have restricted our investigation of flock equilibria to the 
case of translating flocks. In particular, we have not attempted to address 
milling flocks in which individuals travel roughly along circular paths creating 



a vortex-like flock pattern. Some papers investigating this type of flock are 
• flQ 11, 27-29, 29, 31, 38|. 

Finally, we note that field data on starling flocks jlj] show that the density 



Refs 



of birds increases as the border of the flock is approached. This behavior is 
not present in our solutions. Thus our model would need modification to 
describe this feature of starling flocks. Nevertheless, we believe that our 
model is still useful in general, and this may be more strongly the case for 
animal groups of other types than starling flocks, e.g., fish schools. 
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4. Long- Wavelength Instability of Planar Equilibria 

Here we analyze the stability of our planar one-dimensional equilibria ob- 
tained in Sec. |3jto long-wavelength perturbations. For simplicity, we consider 
an equilibrium of individuals with constant flight velocity lying in the surface 
of the sheet, v = v z Q . We now introduce a perturbation depending on x 
and y, 

p(x)=p (x)+5p(x)e sWk y, (41) 
v(x) = v + 5v(x)e st+iky , (42) 

where the density profile po(x) is determined in Sec. 13.41 Using Eqs. ( 14"T|) and 
d32D in Eqs. (pQ) and ©, we obtain 

s5p = -V-(p 5v) (43) 
1 BP 

s 5v = — V (p + Sp) - V (U + SU) - 2t- 1 (v • 5w) v . (44) 

Po op 

For perturbations that are functions of x and y, the perturbed velocity of the 
unstable mode will lie in the x — y plane, and, since we take vo = vqZq, we 
have that the last term on the right-hand side of Eq. (jUJ) is zero (v • 5v — 0) . 
Integrating across the flock, Eq. ([13]) yields 

5pdx = —ik / p 5v y dx. (45) 

■x B J-XB 



The y component of the velocity equation, Eq. (144]) . is 

dx. (46) 



x B r-x B 

Po5v y dx = —ik I 

■x B J -x B 



BP 

po$U + —dp 
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Combining Eqs. f l4"5]) and (j46p we get 



■i'B 



5p dx 



•''B 



dP 

p SU + — 5p 
op 



dx. 



(47) 



For long-wavelengths, (kx-s) <C 1, an expansion in the small parameter (kx-s) 
shows that, to lowest order, at each point y the equilibrium in x (i.e., Eq. f JTTj) ) 
holds. This is essentially equivalent to the statement that for (kx^) <C 1, 
derivatives in x are much larger than derivatives in y and t. However, the 
integral, 

N(y,t)= [ pdx, (48) 



giving the flock number density per unit area transverse to x can depend 
(slowly) on y and t, and the quasi-static equilibrium depends on y and t only 
through N. Thus explicitly considering the equilibrium dependence on N by 
writing p and U as p (x,N) and U (x,N), and setting iV = N + SN, we 
can express the perturbed quantities in Eq. fj47|) as 

(49) 
(50) 



BTI 

SU = e st ~ iky —^5N. 
dN 



Using our definition of the cross-sectional number, Eq. ( HHj) . we obtain from 
Eqs. flSD, (USD, and §D]) 



S *5N = -k 2 



•I'B 



J'B 



Po- 



+ 



5N dx. 



(51) 



9iV dp dN_ 

We re-express the integral J po(dUo/dN) dx as follows. First we start with 
the dimensionless potential equation (see Eq. (TPTj) with normalizations chosen 
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to set £ 2 



u = 1), 



dx 2 



U -U = po- 



(52) 



We take a partial derivative of Eq. ( 152"]) with respect to N, multiply by U , 
and integrate from — xb to xb (recall that the flock density goes to zero at 
the boundaries) giving 



■'■p, 



Uo 



j-ii 



dx 2 



- 1 



dUo 
dN 



dx 



■l-B 



TT 9 P° A 

Uo^dx. 



We also multiply the original potential equation Eq. 
integrate from — xb to xb, to obtain 



(53) 

by OUo/dN and 



dx 2 



1 



Un dx 



ON 



Po dx. 



(54) 



We then subtract Eq. (153]) from Eq. (154"]) . Using integration by parts twice 
cancels all terms on the left hand side. There are no boundary terms in the 
integration by parts since Pq(±xb) = 0. We then have 

dU 



Thus 



■I'D 



Uo dN dX 



dU , 

PO^TTT dX 



dN 



Po dx. 



■J'B 



dN 



l_d_ 

2dN 



J'B 



PqUq dx. 



Using this in Eq. ( j5Tj) gives the final result, 



S 2 SN = -k 2 5N 



where 



_d_ 

dN 



E 



■l-B 



-PqUq + P(po) 



dx = -k 2 5N 



dE 



J'B 



-p Uo + P(p ) 



dx. 



(55) 



(56) 



(57) 



(58) 
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Thus we see that the equilibrium is unstable if BE /ON is negative. For the 
case 7 = 2, we can analytically obtain E and show that dE/dN < by using 
explicit forms for po(x) (Eq. f )25|) ) and Uq(x) (see Sec. 13. 3j) and inserting them 
in Eq. flo^j) . In addition, our numerical calculations of E as a function of N 
for pressures P = cp 1 show that BE / dN < for all 7 where finite equilibria 
exist (7 > 1). See, for example, Fig. [8]for plots of E(N) for 7 = 1.1, 1.4, 1.7, 
3, 5, and 7, as well as P = p 2 for K = 1.25, 2, and 5. These plots show that 
E(N) decreases monotonically with N in all cases. Thus we conclude that 
our planar equilibria are always unstable. 




Figure 8: E versus N for various pressures, a) Numerical evaluation of E versus N for 
P = p 1 for 7 > 2. b) Numerical evaluation of E versus N for 1 < 7 < 2. c) Plot of the 
analytical result of E versus N for P = p 2 . 
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We interpret the basic reason for the long-wavelength instability that 
we have found as being due to the long-range attraction of flock members. 
Since, at long wavelength, attraction cannot be balanced by repulsion, long- 
wavelength density modulations transverse to x tend to grow. Eventually, as 
the density enhancements collapse under the attractive force to smaller spa- 
tial size, we expect that the repulsion will come into play, and an equilibrium 
where attraction and repulsion balance in all directions will be established. 
Our three-dimensional ellipsoidal equilibrium (Sec. 13. 5j) is of this type and is 
thus expected not to be susceptible to this instability. Our planar analysis 
(Sees. I3.1[ 13. 2[ 13. 3[ and 13 .4[) has been useful in gaining understanding and 
in formulating our three-dimensional ellipsoidal equilibria (Sec. 13. 5fl . How- 
ever, based on the above discussion, we expect only our three-dimensional 
equilibria (and not our planar equilibria) to be relevant in modeling real 
situations. 

While the discussion of stability in the previous paragraph suggests that 
the particular mechanism causing long-wavelength instability of the sheet 
equilibrium may be absent in our elliptical equilibrium, we caution that other, 
shorter-scale instability mechanisms are possible. For example, numerical re- 



sults from an agent based model (Ref. 39]), suggest that interactions tending 



to equalize the alignment of flock members may be necessary for stability of 
a finite group. This kind of interaction would not affect our equilibria, which 
have all flock members aligned, but might conceivably be necessary for sta- 
bility of our ellipsoidal equilibrium. This issue deserves further study (e.g., 
one simple way of including alignment interaction in our fluid-like formula- 
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tion is given in Ref. 25j). Another issue that we have not addressed and 
that is open for further study is that of stability in the case of topological 
interactions (Sec. 13. 6j) . 



5. Conclusions 

In this paper we have used a dynamic continuum model to investigate 
animal flocking. We model the short-range repulsion of flock members by 
a pressure-like term, P = cp 1 . Adopting the hypothesis that 'reasonable 
equilibria' go to zero flock density at a well-defined flock boundary, we find 
that such equilibria only occur if 7 > 1, and we investigate the form of 
these equilibria. Adopting the further reasonableness hypothesis that, as the 
spatial size of the flock gets larger and larger, the interior density remains 
bounded, we find that 7 must exceed 2. Furthermore, for 7 > 2 we find that 
large flocks have an approximately constant interior density. 

We considered planar flocks depending only on one Cartesian coordinate 
as well as ellipsoidal shaped flocks. The analysis of planar flocks was analyt- 
ically convenient because the basic nonlinear equilibrium equation could be 
integrated by the method of quadrature. Insight gained from the analysis of 
planar flocks allowed a good understanding of ellipsoidal flocks, particularly 
in the case of large flocks with 7 > 2. 

Our ellipsoidal flocks can be thought of as resulting from anisotropic sens- 
ing and response of flocking individuals to their neighbors. In our modeling 
scheme this resulted in an ellipsoidal shape with principle axes aspect ratios 
that were independent of flock size, and this feature is in apparent agreement 
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with observation ll | . 

Recent work [24j has suggested that an individual in a flock interacts with 
an approximately fixed number of other flock members. Thus the geometri- 
cal interaction range is larger /smaller when the flock density is lower/higher. 



Reference 24( calls this behavior interaction through 'topological distance'. 
We have shown that our model and results based on purely-geometrical 
distance-determined interactions can be simply transformed into a model 
and corresponding results incorporating topological distance. Thus our main 
conclusions are robust in that they apply for geometrical distance and topo- 
logical distance interaction models. 

Finally we investigated the long-wavelength stability of our planar equi- 
libria and found that they were always unstable to perturbations within the 
plane. We interpret this instability as being due to the domination of long- 
range attraction over repulsion at long wavelength. Thus, of our two types 
of equilibria solutions (planar and ellipsoidal), we believe that the ellipsoidal 
type is more relevant for modeling. 

This work was supported by ONR grant N00014-07-1-0734. We dedicate 
this paper to the memory of our coauthor, Parvez Guzdar, who passed away 
recently. 
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